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Abstract 

We show that Peccei-Quinn and lepton number symmetries can be a natural outcome in a 3-3-1 model 
with right-handed neutrinos after imposing a Zn®Z 2 symmetry. This symmetry is suitably accommodated 
in this model when we augmented its spectrum by including merely one singlet scalar field. We work out the 
breaking of the Peccei-Quinn symmetry, yielding the axion, and study the phenomenological consequences. 
The main result of this work is that the solution to the strong CP problem can be implemented in a 
natural way, implying an invisible axion phenomenologically unconstrained, free of domain wall formation 
and constituting a good candidate for the cold dark matter. 
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I. INTRODUCTION 



The Standard Model (SM) of strong and electro-weak interactions, SUc(3) <8> SUl(2) <8> C/y(l), 
has shown its extraordinary accuracy in explaining many features of particle physics along the 
years. Among the issues not covered by this successful model there is the fact that the QCD 
vacuum has a nontrivial structure revealed by its non-perturbative regime, implying the so called 

n 

strong-CP or problem (the subject is widely reviewed in Ref. jy). The violation of CP by strong 
interactions appears in the theory after the introduction of instanton solution to solve the J7a(1) 
problem |2j. It induces the so called #-term in the QCD Lagrangian, which violates P, T and 
CP. Additional electro-weak effects change this term proportionally to Det[M], where M is the 
quark mass matrix. The effective #-term, e ff, is observable through the electric dipole moment 
of neutron, whose experimental bound implies the upper limit \0 e ff\< 10" 9 P. The smallness of 
e ff is what we call the strong-CP problem. 

Among the several solutions proposed to solve the Strong-CP problem, there is one which is 
particularly elegant. It was introduced by Peccei and Quinn and consists of imposing a global 
chiral symmetry, known as Peccei-Quinn (PQ) symmetry, Upq(1), to the classical Lagrangian so 
that the dynamics of the theory sets e ff to zero. Due to the breaking of PQ symmetry a massless 
pseudo-scalar is generated, the axion, which couples linearly to the axial anomaly. When this 
axion develops a vacuum expectation value (VEV), vpq, it produces a further displacement on 
e ff, making it disappear in favor of a dynamical field, the physical axion, eliminating the strong- 
CP violating term of the theory. The breaking of PQ symmetry brings a new scale into the theory, 
fpQ, the axion decay constant, bounded by astrophysical and cosmological data, and its allowed 
range is, 10 9 GeV < f PQ < 10 12 GeV p. 

The first class of models introducing the axion via PQ symmetry in the context of SM were 
obtained by Weinberg and Wilczek (|. This axion was soon shown to be unrealistic mainly due 
to its non-suppressed coupling to light matter fields pj, which happens when vpq is of order of 
the electro-weak scale. Viable models to solve the strong-CP, introducing an invisible axion, were 
devised by Kim and independently by Shifman, Vainshtein and Zakharov Q], the KSVZ axions, 
and by Dine, Fischler and Srednicki as well as Zhitnitskii ^], the DFSZ axions. Both make axions 
invisible by increasing vpq (the larger vpq, the weaker the axion-matter coupling) and obtain 
the axion through a singlet scalar. In the KSVZ axion model the ordinary quarks and leptons 
do not carry PQ charges, some heavy new quarks have to be included which carry this quantum 
number. On the contrary, in the DFSZ ordinary quarks and leptons do carry PQ charges, although 
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these fermions do not couple directly to the singlet, which happens only at the loop level through 
interactions in the potential. 

The possibility of an invisible axion makes the PQ approach even more attractive since in this 
case the axion is a natural candidate for explaining the existence of cold dark matter (CDM) [ic| . 
This is possible because the axion receives a tiny mass through chiral anomaly, rv? a ~ ^qcdI fpo> 
amounting to a mass of 0(1O -5 ) eV. However it is not easy, in general, to obtain the required 
PQ symmetry in a natural way, most models have to impose it from the beginning, weakening 
such a solution to the strong-CP problem. That is the reason we concentrate here in a class of 
models where the symmetry would arise automatical v, namely a particular version of the SUc(3)® 



SU L (3) <g> U Y (1) model (3-3-1 for short) 
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In 3-3-1 models the anomaly cancellation requires a minimal of three families (or a multiple of 
three in larger versions). Besides, there is a bunch of new particles and interactions which make 
these models phenomenologically rich and attractive as an alternative to the SM. If we assume that 
in the realm of intermediate energy there are no exotic leptons, then the 3-3-1 symmetry allows for 
only two possible gauge models for the strong and electroweak interactions, which will be referred 
as version I and version II. 

In the most popular one, version I the triplet of leptons is composed of {vl ,1^ Jr) 1 \ h 
contains exotic quarks with electric charge 4/3 and 5/3, and doubly charged bilepton gauge boson, 
C/ ±=l= , which prompts rare lepton decays. It also implies an upper bound on the Weinberg angle, 
s'm(9y\?) < 1/4. Version II is the 3-3-1 model with right-handed neutrinos |l2j]. In it the triplet of 
leptons is constituted by {vl Jl , u m)J • ^ s bilepton gauge boson is neutral and their exotic quarks 
carry usual charges, 1/3 and 2/3 [la]. 

The physical properties of these models were investigated in several works and their different 
aspects became evident |l^. \\\. Among these differences it is noticeable that version I requires a 
minimal of three triplets and one sextet of scalars in order to generate the masses for all fermions 
and gauge bosons while version II does the same job with only three triplets. 

Their shared aspects include the naturalness of massive neutrinos, with the difference that in 
version I neutrinos are Majorana-type, while in version II they are Dirac-type. Besides, from their 
structure these models dispose of enough constraints upon the U{\)m quantum numbers leading 



to the correct pattern for electric charge quantization [17 1. Another of these aspects is that also 
the PQ symmetry and the leptonic symmetry can emerge naturally in both versions [lsl Ilfll. . 

Since the version II of 3-3-1 was observed to possess the PQ symmetry with a smaller con- 
tent although in that context the axion was of the Weinberg-Wilczek kind, we decided to 
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chose this more economical model and investigate the possibility of obtaining an invisible axion 



by including only one extra scalar singlet fielc 
in version II of 3-3-1 was recently addressed 



in the model. The presence of CDM candidates 



211 ] . but here we wish to have the axion playing 



such a role. There is a crucial issue that has to be addressed when trying to stick with a CDM 
singlet axion though. It concerns the fact that gravitation induces dangerous effective terms in 
the Lagrangian, explicitly breaking any global symmetry of the theory. In particular, focusing on 
Upq(1), this breaking implies a huge contribution to the axion mass. There remains the question 
whether an appropriate mechanism exists in order to avoid such terms, stabilizing the axion. For- 
tunately, the annoying terms can be conveniently suppressed by the presence of suitable discrete 



22 1, that a local con- 



symmetries. Moreover, it was noticed in late eighties by Kraus and Wilczek 
tinuous symmetry at high energies manifests at low energies as discrete (local) symmetries which, 
differently from global ones, are expected to be respected by gravity. This means that the needed 
discrete symmetries can arise in a rather natural way if we assume some underlying local continuous 
symmetry. 

Discrete gauge symmetries have been used to stabilize the axion in a model with extra- 
dimensions by Kamionkowsky et al. 23] more than ten years ago. It was also pointed out that 
large discrete symmetries can naturally arise in the context of string theories |2J|. Also, in an 
attempt to prevent B — L violation in a class of super symmetric standard model, large discrete 
symmetries were imposed, implying an automatic PQ symmetry, stabilized against quantum grav- 
ity effects [25! • For what we are concerned, it was noticed in Ref. ^| that 3-3-1 models possess a 
large enough number of fields to accommodate large discrete symmetries, Z^. And the larger N 
is, the higher are the number of suppressed unwanted terms in the Lagrangian. In order to obtain 
a Z\% symmetry, the authors in Ref. 19] added some extra fermion fields to the model, resulting 
in an automatic PQ symmetry and the axion protected under gravitational mass corrections. This 
constitutes an additional motivation for considering these 3-3-1 models to obtain the invisible axion 
and solve the strong-CP problem. 

This work is divided as follows. We first introduce the model in section [H] In section II I II we 
impose a Z\\ (g> Z<i symmetry, associating the appropriate charges for the fields and obtain that 
the resulting Lagrangian is invariant under Upq(1), identifying the correct PQ charges. This is 
done within the same spirit as that presented in Ref. 19j, |2(| , assigning charges under a discrete 
symmetry group to the fields at hand and observing that a PQ symmetry emerges automatically if 
a Z2 is also imposed. We will see that in this case, also lepton number symmetry arises naturally. 
In section IIVI we analyze the symmetry breaking pattern of the model, recognizing the axion and 
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its couplings. We finally present the conclusions in section IVl 



II. THE MODEL 

Our investigation on this work relies on the version II of the 3-3-1 models Its left-handed 
lepton content comes in the fundamental representation of the SU(3)l, composing the following 
triplet, 



ft 



~ 1,3, 



-1/3) 



(1) 



and the right-handed leptons are singlets, 



e aR ~ (1,1,-1) , 



(2) 



with a = 1,2,3 representing the three known generations. We are indicating the transformation 
under 3-3-1 after the similarity sign, Differently from version I, right-handed neutrinos are 

already present instead of exotic leptons. 

In the quark sector, one generation of left-handed fields comes in the triplet fundamental rep- 
resentation of SU(3)l and the other two compose an anti-triplet with the following content, 



Q 



iL 



-UiL 
\ d 'iL J 



~(3,3,0), Q 3L 



C?3L 
\ U 3Lj 



(3, 3, 1/3), 



(3) 



and the right-handed fields, 



u iR ~ (3, 1, 2/3) , d iR ~ (3, 1, -1/3) , d! m ~ (3, 1, -1/3) 
u 3R ~ (3, 1, 2/3) , d 3R ~ (3, 1, -1/3) , u' 3R ~ (3, 1, 2/3), 



(4) 



where j = 1, 2 represent different generations. The primed quarks are the exotic ones but with the 
usual electric charges. 

In order to generate the masses for the gauge bosons and fermions, the model requires only 
three triplets of scalar s, namely, 



X 



( x o\ 

X~ 

\x'° ) 



~(1, 3, -1/3), rj 



( rf \ 



(1,3, -1/3), p 



\P' + J 



(1,3,2/3). (5) 



With these scalars and matter fields we can write the following Yukawa interactions 

C Y = G^slu'^x + GiiQ iL d' jRX * + Gl a Q 3L u aR r, + G\ a Q lL d aR ri* 

+GfQ 3L d aR p + Gi a Q lL u aRP * + h ab f aL e bRP + K^ifaLUhLf^k + H.c. (6) 

After the breaking of the 3-3-1 symmetry the vector gauge bosons W^, , U° and U ^ interact 



with matter as follows 
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+(u' 3L ^d 3L + u iL ^d! iL )V+ + (u 3L ^u' 3L - d' iL ^d lL )Ul] + H.c. . (7) 

It is through these Lagrangian interactions, C Y and C , that we can recognize particles that 
carry lepton number L such as total lepton number is conserved at this level. Prom these interac- 
tions we have 

L(V + , u> 3 , V '° , p'+) = -2 , L(U° , 4 , x° , X~) = +2. (8) 

Notice that the new quarks, u' 3 and d\ are leptoquarks once they carry lepton and baryon 
numbers; are charged vector bileptons while U° and U ^ are neutral vector bileptons. We have 
also charged scalar bileptons and two neutral scalar bileptons. These last ones would be important 
in studying spontaneous breaking of lepton number, if the associated global symmetry is conserved 
by the potential, leading to the so called majoron, as discussed in Ref. j^ . 

We include also an additional singlet scalar field, <p ~ (1, 1, 0), in order to complete the spectrum, 
allowing for the desired discrete symmetry which will enable us to get an axion protected under 
large gravitational contribution to its mass. 

Finally, we can write the most general, renormalizable and gauge invariant, potential for this 
model. We divide it in two pieces, one hermitean, Vh, and one non-hermitean, Vnh, which can be 
written as, 

V H = pl<p 2 + p\x 2 + p\tf + M> 2 + Aix 4 + W + hp 4 + A 4 (x f x)(^) + A 5 (x f x)(pV) 



+A 6 (r ? t r ? )Go t / o) + X 7 {x ] v){ri ] x) + Mx ] p){p ] x) + p){p ] v) + A 
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+A 11 (<^*)(x t x) + A 12 (#*)(pt p ) + AisCWX^tj) , (9) 

and 

Vnh = p 2 xn X ] r] + hx'v^ + hx^fi* + A i4(xM 2 + Ai5X f ??^ + ^leX 1 V<f>* ^ + Ai7X f #*0* 

+ -/= e ijfc UmPjXk + UViVjPk + hXiXjPk) + e 1 ^ {MmPjXk + ^iWiVjPk + MoXiXjPk) <t> 



+A 26 (x t ??)(P t P) + A 27 (x t ??)(x t x) + H.c. . (10) 

With this at hand we have all the necessary ingredients to associate a discrete symmetry, Zn, 
to the model. This will allow us to eliminate several terms in the non-hermitean potential Eq. 1)10(1 
and verify that we need only an additional Z2 symmetry to have PQ symmetry naturally, assigning 
the appropriate PQ charges to fermions and scalars. 



III. Z u AND PQ SYMMETRIES 

A discrete symmetry Zpj can naturally be accommodate when the theory has enough number of 
fields in its spectrum. It was observed that this is the case for the SM when some scalar multiplets 
and right handed neutrinos are added [2a], or for the minimal 3-3-1 model when only right handed 
neutrinos need to be included It was obtained that a Z13 local symmetry could be imposed 
in this way, leading to a natural PQ symmetry. We remark that this idea was first pursued by 
Lazarides et al. in the context of models embedded in superstring theories j^ . 

Here we are going to apply such idea to the version of 3-3-1 model presented in section [Hj 
that has right handed neutrinos in its fundamental representation. It was observed that an axion 
might be a natural outcome when a Z2 symmetry was imposed in this model [liJ ]. Although this 
axion is of the Weinberg- Wilczek kind, thus phenomenologically discarded [3], if we consider the 
enlarged spectrum with a singlet scalar, eft, the axion can be a mixing of this field with other 
scalars in the model, with its major component being the pseudo-scalar part of the <j) field. Then, 
a discrete symmetry can be imposed allowing for an axion also protected under gravitational mass 
corrections. 

To proceed in this way we first assign the Zn charges to all independent fields, and check for 
additional symmetries appearing after eliminating forbidden terms under Zn . It will turn out that 
a chiral U(l) symmetry arise, and we will see that it is possible to identify it with PQ symmetry. 
It would be interesting to have a Z13 symmetry so as to obtain a PQ scale in its upper limit, 
vpq ~ 10 12 GeV. Although the model disposes of 14 independent multiplets, it is not possible to 
accommodate a symmetry greater than Z\2, because the Yukawa interactions in Eq. © imply some 
constraints over the allowed Zn charges. It is clear that N = 12 is the value of N that allows for 
a maximal protection of the axion under gravitational effects in this model. However, besides the 
seemingly difficulty of avoiding to repeat the phases of the multiplets, the singlet <p would have to 
acquire a very specific phase since twelve is not a prime number. In other words, any even phase 



would make the transformation to belong to a smaller discrete symmetry, jeopardizing our intent 
of suppressing some high order operators involving products. For this reason the largest discrete 
symmetry we can use is Zn, which allows any phase to eft, except the trivial one. 

The effective operators responsible for the gravitational mass contribution are of the form 
(j) n /MpJ A . A Zjy symmetry automatically suppress terms of this kind till some n = N — 1. The 
main surviving term contributing to the axion mass is the one with n = N . It is true that with Z\\ 
the axion is protected only for energy scales not bigger than (<p) ~ 10 10 GeV. Nevertheless, this is 
not a threat for the model since we still have values for the 9 angle and axion mass (gravitationally 
induced) |3(l |. 



V m pi 

e ~ W_ ~ in- 19 fin 

1V1 PI QCD 

where we have used Mpi ~ 10 19 GeV and Aqcd — 300 MeV, and m a ~ 10 -5 eV is the instanton 
induced axion mass. These values are consistent with astrophysical and experimental bounds (see 
PDG |3j). If we had taken (cp) ~ 10 11 GeV, the axion would still be protected under gravitation, 
but the 6 value would be on the threshold of its bound 9 e ff < 10 -9 . So we can have a valid solution 
to the strong-CP problem for Z\\ for scales (0) < 10 10 GeV in this version of 3-3-1. In order to 
seek for this solution let us proceed further by first assigning the correct Z\\ charges to the fields. 
Defining u>k = e m ~ , {k = 0, ±1, ±5} the Z\\ transformations are given by: 



d> - 


-> Ui<f>, 


faL 


-1 t 
^1 JaL , 


p - 


"> U 2 p, 


d a R 




x - 








QiL 




d' iR - 


-1 7/ 


v - 




U a R 






-* ^0 QsL ■ 







(12) 

At this point it is possible to go back to the potential, Eq. ()10|) . and note that this symmetry 
eliminates all non-hermitean terms except three, namely, V'P* 'P* ^ VPX'P-, WP4>*- 
If, besides the Z\\ symmetry we impose a Z 2 symmetry that acts as, 

(<P,X,d' R ,u' m ) -» -(«/», x,d' R ,u' m ), (13) 
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with the remaining fields transforming trivially, the only term which remains in the non-hermitean 
potential is the r/px4>- It should be noted that the Yukawa interactions in Eq. © do not allow for 
terms which interchange x *~~* ?7> since they do not respect Z\\ <g> Z2 given by Eqs. ()12|) and (|13|) . 

We have the stage settled to see that an automatic PQ symmetry arise in the model. To 
achieve this conclusion we start by assigning the PQ quantum numbers such that quarks of opposite 
chiralities have opposite charges, yielding chiral quarks under Upq(1) transformation, 



UaL ~ 


-> e 


~ iaXu u aL , 


U a R ~ 


-> e iaX -u aR 


U 3L ~ 


-> e 




U 3R ~ 


-> e iaK u' 3R 


daL ~ 


-> e 


~ iaXd d aL , 


d a R - 


- e iaXd d aR: 


d'iL - 


-> e 


- iaX '<L , 


d' iR - 


* e^d' iR . 



(14) 



For the leptons we can define their PQ charges by, 



e aL -> e iaX ^e aL , e aR -> e iaX ^e aR , 

v aL - e iaX ^ aL , ^ - e ioX ^ aR . (15) 

With these assignments and taking the Yukawa interactions in Eq. ((BJ) into account, as well as 
the non-hermitean terms rjpx<f>, we easily see that the PQ charges for the scalars are constrained 
and imply the following relations: 

Xd = —X u , Xd> = —X u i , X u = X eR , X e = X vR . (16) 

We can make the further choice X& = X$ , leading to 

Xd = X^ = —X u = —X u i = —X e = X eR = X v = — X vR , (17) 

implying that the PQ symmetry is chiral for the leptons too, and the scalars transform as, 



- 


-> e~ 2iaXd <j), 




^ e 2iaX dr j) 


rf - 


■* V~ , 


r/° - 


+ e 2iaX d ^lO 


p + - 


- P + , 


^° - 


^ & -2iaX d 


p' + - 


-> P + , 


x° - 




x~ - 


-> X _ ) 


x'° - 





(18) 
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Now it is transparent that the whole Lagrangian of the model is Upq{\) invariant and the 
strong-CP problem can be solved in the context of this model. The strong-CP violation angle is 
given by the sum over the quarks PQ charges, which translates to 

6 — ► 6 ± 2aX d . (19) 

This result is possible in this version of 3-3-1 because the PQ charges of the exotic quarks, d! and u' 3 
do not cancel exactly for the case of interest here, X u = —X&- Moreover, the model is particularly 
attractive in the sense it does not present the domain wall problem This means that there is 
no discrete subset of PQ symmetry that leaves the axion potential invariant, i.e., $Zn C Upq(1) 
such that V ax i on {9) is invariant. This is similar to what happens in the 3-3-1 version discussed in 
Ref. [l9(, although there right-handed neutrinos had to be added to the model besides the singlet 
scalar. 

It is remarkable that under Z\\ <g) Z-i, not only the PQ chiral symmetry is automatic but the 
lepton number symmetry also appears naturally in the model, once the possibly non-conserving 
lepton number terms present in the potential completely disappeared. In this sense, discrete 
symmetries originating at some high energy scale seems to be enough to generate the desired 
global symmetries we need at lower energies. 

We finally write the most general potential invariant under 3-3-1 and Z\\ (g) Zi (or Upq{1) and 
Lepton number) symmetries , 

V(ri,p,x) = VH + X^mPjXkt + H.c, (20) 

where Vjj is given in Eq. l|9*|). 

In the next section we are going to use this potential to recognize the axion, the Goldstone 
boson originating from the breaking of the PQ symmetry, and verify that it is constituted mostly 
of the singlet cf). 

IV. SPONTANEOUSLY BROKEN PQ SYMMETRY 

The potential given in the previous section, Eq. (|20|) . allows us to obtain the mass eigenstates 
for the scalars, so we can identify the Goldstones which are absorbed by the massive gauge bosons 
and extract the axion in terms of the interaction eigenstates. To accomplish this, let us consider 
that only x'°j /A V° an d 4> develop a vacuum expectation value (VEV) and expand such fields 
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around their VEV's in the standard way, 

1 



X 



1 



^_ ,v x , + R x , + il x , ) , rj = -j= (v v + R n + iI v ), 



1 



-j={v p + Rp + ilp) , 



1 



(21) 



With such expansion, the next step is to get the constraints that lead to the minimum of the 
potential, 



M x + Any + —v v + — v p + — ^ + -2- = 0, 



2 

Ae 
2 

Ae 
2 



2 , \ 2 , A4 2 , A6 2 , A13 2 , 
Mr, + A2«„ + — U x / + — V p + — U + -2=0, 



■• , v 2 , A5 2 , A6 2 , A12 2 , ^ n 

M P + A 3 v p + —V + — v v + — ^ + ^ = 0, 



2 

Vlj 

y 

Y 



A 



2 1 \ 2 . An 2 , A12 2 , A13 2 , A 
H<t> + M0V4, + —v x , + — v p + —v v + 3 = 0, 



(22) 



2 x 2 p 2 

where we have defined ^4 = X^v^VpVyv^. Substituting the expansion in Eq. 1)211) in the potential 
Eq. (|2U|) and using the constraints above, we get the following mass matrix, M^(R X , R^), for the 
real scalars in the basis, {R x , R^i), 



4 2t^7 4 2v x /v v 

\ 4 2v x ,v v 



* 7 V A 



(23) 



4 ^ 

and the mass matrix, M^(R X > ,R n ,R p , R^), in the basis (R x ' , R v , R p , 
( A x tA 



+ 



o + 







+ 






+ 




A 


2~t?7 

X 


2 




2 


2jj p u x > 


2t)0U x / 


A 


A 2 ^ 




A 




+ 


A 


A 


2d^u x / 






2 


2jj p u,, 








+ 


A 


As^-p 




A 


A 


2v p v x i 


2 


2v p v v 




2^ 


2v p v^ 



\ 



2v v d _ 



(24) 



2v 4> v x , 2v n v (f> —p-cp t ^ 

These bases are not coupled, that is the reason we have two squared mass matrices. From the 
first matrix, Eq. (|23|) . after diagonalization it is easy to recognize the following massless scalar in 
it, 

1 



Rg 



vjj + v 2 x , 



(v n Rjf - v x >R x ) 



(25) 



The other real scalar mass eigenstate is orthogonal to this one and those coming from the diag- 
onalization of matrix Eq. (|24|). which are a little more intricate but fortunately we do not need 
them for our purpose. 
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Regarding the pseudo-scalars, similarly to the real scalars, we obtain the following mass matrix, 



Mj ( I x , I„> ) in the basis (I x , L, 



rf ) ) 



V 



A 7^ 



\tv x iv v 



A 

x 

+ 



A ~T 



X 



A 



(26) 



and the mass matrix, Mj(I x > ,I n ,I p , 1$), in the basis (I x > ,I v ,I p , I^), 



( i 



.4 
1 



i 



i 



v v v x , v p v x i v x iv$ 



1 



1 



1 



v P v x i V v Vp 
1 1 



1 



1 



1 

9 



(27) 



/ 



From these matrices we can easily obtain the Goldstone bosons and identify the axion as the one 
whose main component is in the <f> direction. The Goldstones and the pseudo-Goldstones are listed 
below, 

1 



Gi 
G 2 

G 3 

PS! 

PS 2 




v<t,v x , 



V P V 1 V X' 



v P vy x , 



v p v 4> v x > 



1 + 
i 



Vj. v\ 
+ 7V + t£- 



x' 



l<j> H i x H ir) H 1 p 



2 



7J' 



(28) 



where we have defined £ = fp(f^ + u^/). In the above equation, the axion is identified as a, and 
the remaining three Goldstones, G±, G 2 and G3, together with the Goldstone in Eq. (|25|) are 
those eaten by the 4 neutral gauge bosons of the model. The last two linear combinations of the 
interaction states in Eq. (|28jl. PS\ and PS 2 , are the massive eigenstates or pseudo-Goldstones. The 
important point that can be extract from these results is that our axion has a small component of 
I x '. Since, ~ 10 10 GeV and ~ 10 3 GeV, this component is very suppressed and, as expected, 
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our axion is invisible being almost exclusively the imaginary part of the singlet (p. Besides, since 
I x i couples only to the exotic quarks, our axion is very different even from that obtained in the 
version of 3-3-1 in Ref. ^{J, which does couple to neutrinos at tree level. Its coupling can be easily 
obtained after rotating the mass eigenstates, Eq. (|28|). in terms of the interaction eigenstates, and 
it translates into the following Lagrangian term, 



£aq'q' 



Gtu'^-G^^dU a + H.c, (29) 



which are very weak for G\,G % 2 3 ~ of the order of unit, since v x >. The pseudo-Goldstones, 

PS\ and PS2, are more strongly coupled to fermions and, differently from the axion, also couple 
to ordinary matter. This leads us to conclude that, in this model, the only candidate for cold dark 
matter is the axion. We could check if the real massive scalars could fit for this role also, but 
a rough numerical approximation just confirms that they behave as their partners, as we could 
expect. 

We also checked the coupling of our axion with photons. It is defined through the effective 
Lagrangian term, 

C^ = ^L^F^. (30) 

327T J VpQ 

In the present model only exotic quarks participate in the loop leading to the above anomaly term, 
which leads to 

Ca 77 = E X <Q* = H ~ - 44 • ^ 

Sv PQ j 9 

This value is very similar to those obtained in different models present in literature and can be 
used to make the relevant computations involving axion in astrophysical processes. 



V. SUMMARY AND CONCLUSIONS 



We studied the consequences of discrete symmetries in the version of 3-3-1 model with right- 
handed neutrinos. One of the main points in this work is that global symmetries appear auto- 
matically as consequence of such discrete symmetries in this model. It turned out that, when the 
model has a Z\\ <g) Z2 symmetry, the whole Lagrangian is invariant under U(l) transformations 
and also total lepton number is conserved at the classical level. It is remarkable that this happens 
in this more economical version of the model by adding only one singlet scalar, no other fields are 
necessary, which makes it a suitable model for implementing the strong-CP problem solution. 
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We then recognized the global symmetry identifying it with a chiral PQ symmetry, Upq(1). 
In general, solutions to the strong-CP problem through PQ mechanism lead to the domain wall 
formation, which is a threatening feature to Cosmology, but fortunately model dependent. In this 
version of 3-3-1 model this problem is absent due to the fact that we chose a relation among PQ 
charges, namely Xj = X' d , which avoid this situation. Nevertheless, we have to remark that other 
possibilities would be allowed if we had not imposed a Z2 symmetry. In this case, we could have let 
X' d free and, working with the non-hermitean potential terms, look for the relations among the PQ 
charges that would keep the Upq(1) invariance. Such a relation exists and is given by 3X' d = —X^ 
which is consistent with PQ invariance for all three Z\\ invariant non-hermitean terms, x v4>* ft* \ 
rjpx(f> and rjr}p(j)* . In this situation we would have to address the formation of domain walls, and 
for this reason our previous choice seemed physically more appealing. 

We proceeded with the spontaneous breaking of Upq(1) and studied its consequences, obtaining 
the axion and showing that it is mainly constituted of the singlet. It was shown that it interacts 
with exotic quarks only, through a very suppressed coupling. Hence, our axion is an invisible 
one. However, since global symmetries are not stable against gravity effects, our axion could be 
in danger, receiving large mass corrections and losing its appealing as CDM candidate. This was 
circumvented by assuming that Z\\ is a subgroup of an underlying gauge group at some high energy 
scale. Annoying gravity induced terms contributing to the axion mass, are conveniently suppressed 
thanks to the local Z%\, stabilizing the axion. Although, larger discrete symmetries would lead 
to a better stabilization of the axion against gravity, as well as a 9 e ff safer from experimental 
constraints, we saw that is still possible to have a Z\\ leading to small mass corrections and 8 e ff 
below the bound eff < 1CT 9 if = {(f>) ~ 10 10 GeV. Besides, the fact that 11 is a prime number 
allows 4> to acquire any charge under Z\\ except the trivial one. This implies no need of assigning 
a very specific charge to <j) in order to avoid restriction to smaller subsets of Z^ which would not 
lead to axion stabilization. 

Finally, there is a point that has to be highlighted when considering the 3-3-1 version with 
right handed neutrinos. In the form it was presented here the model generates arbitrary masses 
for two neutrinos only, which can be deduced from the anti-symmetry of the Yukawa coupling h a i,/ 
in Eq. Q. This leads to a anti-symmetric neutrino mass matrix, implying a zero eigenvalue and 
two degenerate ones. Although the massiveness of neutrinos is not the issue here, it would be nice 
to have a model that at least produces an appropriate arbitrary mass pattern for all fermions of 
the theory. To accomplish this we have to devise some way of eliminating such mass degeneracy. 
We can suggest two ways of doing that. One of them makes no enlargement of the spectrum and 
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seems the preferable one, dealing only with the vacuum of the theory. In Ref. 20], the breaking of 
leptonic symmetry could be achieved only through a non-conserving PQ symmetry term, namely 
the ijpx, when rj' acquires a non-vanishing VEV. However, we have seen that our approach allows 
for an equivalent term which could lead to lepton number violation too, which is r]px<p, with the 
difference that in this case PQ symmetry is still conserved. This term allows for a Majorana 
neutrino mass through loop corrections, making possible a mass matrix which is arbitrary enough 
to accommodate non-degenerate and non-zero neutrino masses. Another way out could be traced 
by including a singlet neutrino in the spectrum. Such neutrino carries the exact quantum numbers 
to provide the required invariance under Z\\ ® Z<i and Upq(1) and leave only non-degenerate 
massive neutrinos in the theory. This second possibility sounds appealing since not many fields 
can be introduced without jeopardizing the desired discrete and global symmetries here studied. 
Whatever nature's choice, both would be adequate to fit in our approach. 

Summarizing, we obtained automatic Upq(1) and lepton number symmetries by imposing a local 
discrete symmetry in an economic version of 3-3-1 model with right-handed neutrinos. We got an 
invisible axion stabilized under gravitational mass corrections, absent of domain wall problem, 
solving the strong-CP problem and constituting a strong candidate to CDM. This results are a 
remarkable achievement of our work, considering that the only additional ingredient we have used 
was the inclusion of a singlet scalar in the model. 
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